-
Notifications
You must be signed in to change notification settings - Fork 26
Fitting new lattice-constrained NP model, PDFs from correction files through helicities #640
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
base: main
Are you sure you want to change the base?
Conversation
| ) | ||
| pdf_hist += "UncertByHelicity" | ||
| pdf_corr_hist += "UncertByHelicity" | ||
| symmetrize = "quadratic" |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Where is this variable used?
| cols, | ||
| **info, | ||
| ) | ||
| if theory_helpers.get("pdf_from_corr") is not None: |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Maybe better:
for pdf_name, pdf in theory_helpers.get("pdf_from_corr", {}).itmes():
logger.debug(
f"Make PDF (from correction file) uncertainty by helicity histograms for {dataset_name} and PDF from correction {pdf_name}"
)
add_pdfUncertByHelicity_hist(
results,
df,
pdf,
pdf_name,
pdf_name,
axes,
cols,
**info,
)
| # alphaS from correction histograms only available for these sets, | ||
| # so fall back to CT18Z for other sets | ||
| if not ("MSHT20" in asname or "CT18Z" in asname or "MSHT20an3lo" in asname): | ||
| if not ( |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
can we dynamically check if the correction histogram is there and otherwise use the default? E.g.
if not asname in theory_corrections.valid_theory_corrections()
| pdf_hist = pdfName | ||
| pdf_corr_hist = ( | ||
| f"scetlib_dyturbo{pdf.upper().replace('AN3LO', 'an3lo')}VarsCorr" | ||
| if self.corr_hist_name == "scetlib_dyturboCorr" |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Should this not be something with "MSHT20" here?
related https://gitlab.cern.ch/cms-wmass/wremnants-data/-/merge_requests/90